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Abstract 

The Martin-Siggia-Rose (MSR) functional integral technique is applied to the 
dynamics of a D - dimensional manifold in a melt of similar manifolds. The 
integration over the collective variables of the melt can be simply implemented 
in the framework of the dynamical random phase approximation (RPA) . The 
resulting effective action functional of the test manifold is treated by making 
use of the selfconsistent Hartree approximation. As an outcome the general- 
ized Rouse equation (GRE) of the test manifold is derived and its static and 
dynamic properties are studied. It was found that the static upper critical 
dimension, d nc = 2D/(2 — D), discriminates between Gaussian (or screened) 
and non-Gaussian regimes, whereas its dynamical counterpart, d uc = 2d nc , 
distinguishes between the simple Rouse and the renormalized Rouse behav- 
ior. We have argued that the Rouse mode correlation function has a stretched 
exponential form. The subdiffusional exponents for this regime are calculated 
explicitly. The special case of linear chains, D = 1, shows good agreement 
with Monte-Carlo-simulations (MC). 
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I. INTRODUCTION 



There has been recent interest in the dynamical behavior of polymers (or more generally 
polymeric manifolds) in a quenched disordered random medium |T]j2|] and in a melt |3|-Q. 
The dynamics of flux-lines in a type II superconductor || also belongs to the same class 
of problems as in [T],[J . The derivation of the equations of motion for the time correlation 
functions has been carried out by making use either of the projection formalism and the 
mode-coupling approximation P-I7J or the Martin-Siggia-Rose (MSR) functional technique 
and the self consistent Hartree approximation |EB|||. The latter theoretical approach (which 
is equivalent to the Hartree type approximation) was earlier also successfully applied to the 
investigation of static properties of different models with |T0|-|l2l or without Hl3Hl5[l replica 
symmetry breaking. 

The basic description for polymer dynamics in general is the so-called Rouse model 
T6|jr7|1 , where the polymer configuration is expressed in dynamical modes. The physical 



background of the Rouse model is very simple: It corresponds to a non interacting chain 
stirred by a white noise random force, in the usual Langevin sense. It is well known from 



experiments JT8[, but also very surprising that this Rouse model provides a good description 
for the melt of the relatively short chains N < N e . At higher degrees of polymerization 
the dynamics of the long chain can be described by the reptation model [|T^,[T^,|rjJ . For the 
case of the short chain melt, it is not obvious that the collisions of the surrounding chains 
close to a test chain add up to a white noise force. On the other hand, the obvious thought 
to explain the Rouseian behavior in short chain polymer melts is first the excluded volume 
screening and second the inactivity of topological constraints ("entanglements") at these 
length scales. The chain in a melt has Gaussian statistics, i.e., R oc \fN can be explained 
by the screening of the excluded volume interactions [0,0 . In this case the chains (or 
generally speaking polymeric manifolds) strongly interpenetrate. The screening of exclude 
volume reduces the upper critical dimension of the interactions (from four to approximately 
two in the case of linear chains), and we can expect that for higher connected polymers such 
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as D - dimensional manifolds the screening becomes also dependent upon the embedding 
space dimension d and the connectivity. 

Nevertheless we show below in this paper that in the case of dynamics the situation is 
more complicated, even for linear chains. The screening of the excluded volume interactions 
leads by no means automatically to the Rouse dynamics even for the short chains. Indeed 
we will show that the interactions introduce a new dynamical regime in 2 < d < 4. This new 
regime is derived on different grounds than those proposed by Schweizer f||| . Moreover we 
need to resolve the following questions: How does the bare monomeric friction coefficient 
£o renormalize due to the interactions of the test chain with the surrounding melt? Under 
which conditions is renormalization relevant? A short presentation of the results of this 
paper was given in a recent letter [^1[ . 



We use here the MSR-functional technique and the selfconsistent Hartree approxima- 
tion for the investigation of the static and mainly dynamic properties of a polymeric D - 
dimensional manifold |22|] (or a fractal) in the melt of the same species. One of the main 
results of the present study is the derivation of a generalized Rouse equation (GRE) for 
such D-dimensional manifold. In this equation the static and dynamic parts are treated 
on a equal footing and in the static limit the screening and saturation of D - dimensional 



manifolds are reproduced in a different way than in p3 . 

We should stress that the manifolds in our consideration have phantom springs, i.e. are 
crossable, so that entanglements cannot occur and the reptational dynamics is not considered 
here. The reptational dynamics is driven by topological constraints and will be considered 
in a subsequent publication. We describe below the manifolds only in terms of connectivity 
and excluded volume interactions. The connectivity defines the D - dimensional subspace 
which is embedded in the Euclidean space of d dimensions. 

This model includes the cases between linear polymer chains, which correspond to D = 
1, and tethered membranes (D=2). By analytic continuation to rational numbers of the 
spectral dimension statements on polymeric fractals can be made. Branched polymers (and 
percolation clusters) correspond closely to the spectral dimension D — 4/3. In a series of 



papers [£3J one of us has considered the different regimes in static scaling. Here we will 
find that these regimes besides the Gaussian one are unstable. Nevertheless we will show 
below that a new dynamical regime for the motion of manifold segments appears. The 
whole dynamical consideration results in a sub-diffusive behavior and exponents, which are 
confirmed for a melt of linear chains in 3 - dimensions by Monte-Carlo (MC) numerical 
investigations [p4H 28f (see also f29 H3lf1 ). 

The paper is organized as follows. In Sec. II, we define the system and derive the GRE. 
To do this we integrate over the matrix collective variables. The resulting action in terms of 
the test manifold variables will be treated in the framework of the Hartree approximation. 
In Sec. Ill, on the basis of the GRE the static and dynamic properties of the test manifold in 
the melt are investigated systematically. We discuss the concept of the static and dynamic 
upper critical dimensions and calculate explicitly the dynamical exponents. Discussions and 
conclusions can be found in Sec. IV. In the light of the GRE some criticism of the Polymer 
Mode Coupling Approximation (PMCA) by Schweizer ||-[7| is relegated to Appendix A. 



II. DERIVATION OF THE GRE 

A. Integration over the collective variables 

We consider a melt of D - dimensional manifolds which is embedded in the d - dimensional 
space. The test manifold is represented by the d - dimensional vector R(x, t) with the D - 
dimensional vector x of the internal coordinates which labels the beads. The total number of 
beads is given byAf = NxNx...xN(D times) . In the same way the manifolds which belong 
to the surrounding matrix are characterized by r^(x, t) (p — 1, 2, . . . , M). The notations are 
taken in such a way, that the boldfaced characters describe the external degrees of freedom 
in the d - dimensional Euclidean space, whereas the arrow hatted vectors correspond to 
the internal D - dimensional space. The model of the melt of M (monodisperse) tethered 
manifolds used in the following is based on the generalized Edwards Hamiltonian 
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H = ^£ dPx (V^^(x)) 2 + t J2 d D x d D x'V 



p=l " " p,p'=X ' 

where e = dT/l 2 is the elastic modulus with the Kuhn segment length I and we consider 
units defined such that the Boltzmann constant ks — 1. 

An additional test manifold is immersed in this melt and is described by the variables 
R(x, t) . As a result the Langevin equations in Cartesian components j for the whole system 
have the form 
d 

^—Rj(x, t) - eA x Rj(x, t) 



+ 



+ 



S 



J d D x'V(R(x,t) -R(f ,t)) 

M . 

<SBi{ V f) Ej d D x'v(R(x,t)-r^(x>,t)) = f 3 (x,t) 



5Rj (x, t) 
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(2) 



and 



£o^rf (x,t)-eA x rf (x,t) 

X M 



+ 



+ 



5rf(x,t)£ 
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Srf\x,t) 



J2 f d D x'V (r (p) (£, t) - r< m >(f , t 

■71=1 

J d D x'V (r (p \x,t) - R{x',t)) = fj(s,t) 



(3) 



where £o is the bare friction coefficient, V(- • • ) is the excluded volume interaction function, 
Ag denotes a, D - dimensional Laplacian in internal space and the random forces have the 
standard Gaussian distribution. 

We find it more convenient to reformulate the Langevin problem (0)-(|3]) in the MSR- 
functional integral representation ||. The generating functional (GF) of this problem can 
be written as 



Z{--} = J DR^x^DR^x.t) 

x exp | — H Rj (x, t) , Rj (x, t) — A\ R(x, t) , R(x, t) } 



(4) 



where dots imply some source fields and the action functional of the test manifolds is given 
by 



A 1 R(f,t),R(x,t) = - / d u x / dtliRj(x,t) 



- / d u x'iR j (x,t) 



T£ iRj{x,t) 
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5Rj(x, t) 



d 

£ogj;Rj(x, t) - ek x Rj(x, t) 



V [R(x, t) - R(f, t)] 



(5) 



Here the hatted vectors R describe the auxilliary (response) fields in the standard MSR 
technique. Their correlation functions are given by < RjRi >= (Causality) and % < RjRi > 
is a response function. 

The influence functional in eq. @ has the form 

. M 

R,R] = - In / l[Dr {p) {x,t)Dr {p \x,t) 
J P =i 



x exp < - A 2 



r (p) 



+ X] / t^ardV / dtiRj(x,t) V (R(z,*) - r^(f',t)) 



p=l 

M 



+ E/ d D xd D x' J dtiff(x,t) V (r^\x',t) - R(x,t 



(6) 



p=i- ' 5r) P) (x',t) 

where the summation over repeated Cartesian component indices is implied and the matrix 
manifolds action is defined 



M 



A, 



r ip \x,t),r (p \x,t)\ = -]T / d D x j dt\irf'(x,t) 

P =i 

M 



to^ P \x,t)-eA x r ( f\x,t) 



m=l 



x: / d D x'iff\x,t) 



5 



5r 7 (p) (f,t) 



-V 



r (p \x,t) -r (m) {x',t) 



+ T£ ir^\x,t) 



(7) 



The representation @-(0) is very useful for performing transformations to collective vari- 
ables as well as integration over a subset of variables. In our case we make the transformation 
to the matrix density p(r,t) and the matrix response field density 7r(r, t): 

M 



p(r, t) = jr j d D x5(r - r ip) {x, t 



(8) 



7P 



M d . 

r >*) = E E / ^ ^(f, *)V i( 5 (r - r (p) (f, t)) . 



(9) 



P=ij=i ' 

The transformation to the collective variables has already been used before for the dynamics 
of semi-dilute polymer solutions |32] and melts [|S|,|3~3f . The principal aim now is to integrate 
in the influence functional (|6]) over the collective variables and as a result to find 

the effective dynamic of the test chain/manifold in the matrix. 

The representation of the influence functional in terms of the collective fields p and it 
has the form 



R, R 



In / Dp(k,t)Dn(k,t) 



xexplw\p,*]-Jdt J ^7r(-M)p(M)V(k) 



+ / d D x I dtiRj(x,t) 



d d k 



d D x I dt 



d d k 



ikjV(k)p(k, t) exp{zkR(£, t)} 



(2tt)« 



vr(k, t)V(-k) exp{ikR(f , t)} 



with a potential W that depends only on properties of the free system 

. M 

W{p, 71} = In / [] Dr^Dr(p) exp {A [r (p) , r (p) ] } 
p=i 

x 5 L(r, t) - J2 J d ° x 5 ( r - r(p) ( f > *)) j 

( M d \ 

x 5 7r(r, t) - E E / *) V i 5 ( r - r(p) (^ *)) 



and A the free system action 
^4o 



r (p) f-(p) 



M 

P =i J 



t ® (P) _ A (P) 



/r 



(p) 



(10) 



(11) 



(12) 



The Eqs. (|TT|)-(|l2j) can be brought into a very compact form in terms of the 2 - dimensional 
column variables p a and Xa 



p a (k,t) 



vr(k,t) 
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(13) 



and 



Xa(M) 



/ d x iRj(x, t)ikjV(-k) exp [— ikR(x, t)] 



V 



/ d D x V(k) exp [-ikR(x, t)} 



(14) 



where a = 0, 1. 

By making use of eqs. (HD-flOD in eq- © we arrive at 



Ft, R, 



In / JjDp a (k,t)exp^{p a } 



r / dt 



d d k 
(2tt) c 



p a (k, t)i7 Q/3 (k)p /3 (-k, t) 



+ 



/^/^X a (-k,t)p a (k,t)| (15) 



with the 2x2- interaction matrix 



1 V(k)\ 



(16) 



\^(k) j 

In eq. (15) the summation over repeated Greek indices is implied. 

Up to now all calculations have been exact. In order to proceed we should specify the 
"potential" W {p a }. Unfortunately the exact form for W is not known explicitly, but in 
ref. [0 the systematic expansion in the p a - variables was given for the first time. Here 
we will use this expansion only up to the second order, which corresponds to the dynamical 
RPA 



W { Pa } = W{( Pa } } + I J dtdt'-0^p a (k, t)W®(k; t - t')p p (-k, ?) + ... (17) 



where 



w$(M) 



(18) 



and (F r ) a/3 (k, t) stands for inversion of the matrix of the free system 



F aP (M) 



f F 00 (k,t) F 01 (k,t)\ 



F 10 (k,t) 
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(19) 



/ 



In the 2x2- matrix fllTf ) F 0Q (k,t), F i(k, t) and F w (k,t) are time correlation, retarded 
response and advanced response functions correspondingly. 

The Gaussian approximation in eq. ([D]) corresponds to the dynamical random phase 
approximation (RPA) [Q. The RPA makes the integration over p a in eq. (15) analytically 
amenable and as a result for the GF we have 



Z{...} = J DR(x, t)DR(x, t) exp {-A [r, r] } 



(20) 



where 



.4 



R, R 



Aq r, r 



- - / dldl'iRjd) 



d d k 



(27T) 



jkjk p \V(k)\ 2 exv{ik [R(l)-R(l')]}ii2 P (l')S'oo(k,t-0 



+ | dldl' i^-(l) y -^^.^(k)! 2 exp {/k [R( I I R( 1')] } -s;,,.( k. / /') 



dldl'iR^l) 



d d k 
(2tt) c 



i^y(k) exp {ik [R(l) - R(l')]} 5(t - t') (21) 



where we used a short hand notation in the form 1 = (x, t), and Aq is the action of the free 
test manifold 



A 



R, R 



dln-RjOL) 



(22) 



5*oo(k, t) and Soi(k, £) are the corresponding elements of the dynamic matrix in RPA, which 
is given by 



Saf,Qt,t)= [C7+(F- 1 )]J(k,0- 



(23) 



The GF ( pO"D determines the dynamics of the self-avoiding test manifold, which is modulated 
by the melt fluctuations given in the dynamical RPA.Q 

To go beyond RPA anharmonicities in the expansion ( fL7|) must be taken into account 



and follow the renormalized perturbation theory which was worked out in ref. [34] for the 



1 In ref. Q the self-avoidance was dropped and instead of RPA-correlators the full correlators were 
introduced. 
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glass transition problem. In case the temperature is much higher than the glass transition 
temperature the fluctuations in a polymer melt are described by the RPA reasonably well 



B. The selfconsistent Hartree approximation 

The resulting action (21) includes the test manifold variables in a highly non-linear way. 
In order to handle this we use the Hartree-type approximation [1,8,9]. In the selfconsistent 
Hartree approximation the real MSR - action is replaced by a Gaussian action in such a 
way that all terms which include more than two fields Rj(x,t) or/and Rj(x,t) are written 
in all possible ways as products of pairs of Rj(x,t) or Rj(x,t), coupled to selfconsistent 
averages of the remaining fields. As a result the Hartree - action is a Gaussian functional 
with coefficients, which could be represented in terms of correlation and response functions. 
The calculations are straightforward and the details can be found in the Appendix B of 
the ref. [8]. The only difference is that here we deal with the self-avoiding D - dimensional 
manifold (see the last term in eq. (21)) in the d - dimensional space and the collective 
dynamics of the melt is treated in the framework of the dynamical RPA. The resulting GF 
takes the form 

Z{---} = J DRDRexp j - A [R,R] 

/poo pt 
dPxdPx 1 \ dt dt' iRj(x,t)Rj(0? \t')X(x 1 x';t,t / ) 
J —oo J—oo 



I —oo J—oo 

COO ft 



/poo pt 
dPxdPx* j dt I dt 1 iRj{x,t)Rj{x,t)\{x, x;t,t') 
J—oo J—oo 

+ - / d D xd D x' / dt dt' iR j (x,t)iR j (x',t')x(x,x';t,t')\ (24) 

2 J J—oo J~oo J 



where 



A(f,f';M') = ~G{x,x' ;t,t') J ^fc 4 |\/(k)| 2 F(k; x, x'; t, t')S 00 (k; t, t') 

- / ^ d k 2 [\V(k)\ 2 S 01 (k;t,t')-V(k)5(t-t')]F(k;x,x';t,t^ (25) 

and 
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X (x,x'-t,t>) = J ^fc 2 |^(k)| 2 F(k;x,x';t,t')5 o(k;t,t')- (26) 
In eqs. (|25| ) -([26| ) the response function 

G(x, x; t, t') = (iR(x, t')R(x, t)) (27) 

and the density correlator 

F(k;x,x';t,t') = exp f; t, t') | (28) 

with 

Q{x, x'; t, t') = (R(f , t)R(x, t)) - (R(f , t)R(f , t')) (29) 

are specified. The pointed brackets denote the selfconsistent averaging with the Hartree-type 
GF (26). 

It is obvious that for the case under consideration the time homogeneity and the 
Fluctuation-Dissipation Theorem (FDT) hold, then 

d 

G(x,x';t-t') =T- 1 —Q{x,x;t-t') t > t' (30) 
S 01 (k;t-t')=T- 1 ^Soo(k;t-t') t > t' . (31) 



By making use of eqs. ( p0|) and ( pT|) in eqs. (24)- (^) and after integration by parts with 
respect to the time argument t', we obtain 
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Z{ - } = J DRj(x, t)DRj(x, t) 



rN roc 

x exp| J d D xd D x'J dtdt' iRj(x,t) 



£ 6(t - t')S(x - x') + 



+ e(t-f)~ / 0- d k 2 \V(k)\ 2 F(k;x,x';t-t')S oo (k;t-t') 



o 



X 



N roo 

d D x d x' / dt iRj (x, t) 



e5(x-a?)A x + J -^k 



\V(k)\ 2 -S st (k)-V(k) 



T 



-F st (k; x, X* ) — 5(x — x 1 ) \ dPx F st (k;x, x 

Jo 



Rj(x, t) 



rN roo 

+ T d D x d D x' / dtdt' 

JO J -co 



f I (0i fc2 |^( k )| 2i71 ( k ^^^-O5oo(k;t-t' 



£ 8{t - t')5{x - x) + Q{t - t') 

iRj(x, t)iRj(x', t') 



(32) 



where -F s t(k; x, x') is the static density correlation function. 

The generalized Rouse equation (GRE), which directly follows from the GF (0), has 
the form 

d „ _ . s r , n , r* ,., . , A d 



£ -R J (x,t)+ J d D x' J dt'T(x,x l ;t-t , )—R j (x',t l 



J d D x'tt(x, x)Rj(x, t) = T(x, t) (33) 



with the memory function 



T(x, a? ; t) 



1 r d d k 



TJ (27T) d 

and the effective static elastic susceptibility 



k 2 \V(k)\ 2 F(k;x,x';t)S 00 (k,t) 



(34) 



fl(x, x 1 ) = eA(x — x')5 x 



d d k 
(2tt) c 



k 2 V{k) 



F st (k; x'x') - 5{x - x') J d D x"F st {k; x, x") 



(35) 



and the random force JF(x, t) has the correlator 



(F(x, t)F(x' , t')) = 2TS ij [£ 6(x - x')8(t - t') + G(t - t')T(x, x'; t - t')} . 



In eq. (35) the effective interaction function 



(36) 
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V(k) = V(k) 



1 - ^V(k)S st (k) 



(37) 



gains the standard screened form [16 



l + V(k)F^>(k)/T 

(where -F s °(k) is the free system correlator), if the standard RPA- result is used for the melts 
static correlator 

Sst (k) = . (39) 

l + V(k)F^\k)/T 

The GRE (33)- (|36|) is the generalization of the corresponding equations given in |§ for 
the case of a test manifold with self-excluded volume. On the other hand the collective 
dynamics of the melt is treated in the framework of the RPA. This is a good starting point 
for a simultaneous consideration of the static and dynamic behavior of the test manifold. 
One should expect, for example, that the reactive and dissipative forces in eq. (33) are 
screened out in different ways. As we will show in the next section this is really the case. 
Explicitly stated, the test manifold could be statically Gaussian, but dynamically it could 
follow a renormalized Rouse dynamics. 



III. STATIC AND DYNAMIC BEHAVIOR OF THE TEST MANIFOLD 

For the following discussion, it is convenient to perform a Fourier-transformation with 
respect to the variable x. We define e.g. the Rouse-mode time correlation function 

C(p,t) = -j j d D xC(x,t)e- l ^P (40) 

and its inverse transformation 

C{x,t) = J d D P C(p,t)e i2 ^ (41) 

where M = N D is the total number of "monomers" (or beads). 
Then eq. (33) leads to the result 
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£o^C(p, t)+Mj Q dt>T(p, t - t')^C(p, + n(p)C(p, t) = (42) 

where T{p, t) and £l(p) are the Fourier transformations of the memory function fl34|) and the 
susceptibility (35) respectively. Below we shall analyse both of them explicitly. 



A. Static properties 



The static limit t — ► + of eq. fl4"2]) can be implemented, if we take into account 



|c(f,») 



T£ Q G(x,t^ 0+) = -dTS{x) 



(43) 



Jt-»o+ 



where the FDT (|30| ) as well as the initial condition for the response function (see eq. (31) in 
H or eq. (3.12) in |l|]) have been used. Then the formal solution for the static Rouse mode 
correlator, C(p) = (R(p)R(— p)), becomes 

d 



C(p) 



(¥) 



(44) 



The eq. ([44]) has the Dyson-like form where the "self-energy" is given by 



V(k)/T 



[F st (k,p)-F st (k,p 



+ V(k)F^(k)/T 

The static correlator C(p) is parametrized by the wandering exponent ( (see eq. 



(45) 



Q s t(x) = I d D p 



1 - e~ l K 



C(p) oc x 2C - . 



(46) 



In its turn the test manifolds static correlator in eq. ([45]) is 



^(k;p) = ^/^exp|-^ 



k 2 l 2 2e 2tt 

, x ^ — « — a; • p 
2d N ' 



(47) 



As a good approximation for the free system correlator, (k) , we will use the Pade formula 

,(0Wn f>M 



FA UJ (k) = .., 

1 J l + iklpAf^D) 

where p is the averaged bead density and 



(48) 
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(49) 



with the Gaussian fractal dimension d°f = 2D/(2—D), the Gaussian exponent Co — (2— D)/2, 
S D = 2tx d ' 2 /T(D/2) the volume of the unit spere and the gamma-function T(x). The system 
of eqs. (|44|) -(f49|) can be analysed self-consistently. 

Let us start from the calculation of the "self-energy" (fH). Because (3VF^ ] > 1, the 
effective screened interaction is proportional to (k). Then using eqs. ([47f )-([49l) in 

eq. ([45D and performing the integration over k first yields 



E(p) = E 1 (p) + E 2 (p) 



(50) 



with 



Si(p) = - Cl (-j 



S 2 (p) 



- 1 n) 



where C\ and c 2 are given by 



S d S D (2d)^ (d 
Cl l 2 (2n) d ' 2(pl d ) \2 



x dtt D-l-«d+2) 



(irp) D Jo 



and 



s d s D (2d)^ 1+i (d+d° f : 



l 2 {2 



It 



2 (pi" 



-p 



x / dtt - 1 -^}^ 
Jo 



(51) 



i-r,f) [}yj^2t) 



(52) 



i-ri§) Q V J^(2*) 



(53) 



Here Jd/2-i(%) is the Bessel function. 

We assume that p = 0(1) but (p/N) <C 1. Physically, the condition for the exponent £ 
comes from the balance between the entropic and the interaction terms in the denominator 
of eq. (fHp. But one should also be wary about the self-consistency condition ptj| ) otherwise 
the result could be different. 
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i. Let us, e.g. assume that 



C(d + 2) < C(d + d°+2)-D. 



(54) 



In this case Si can compensate the entropic term proportional to (np/N) 2 in eq. 
and one should claim ((d + 2) = 2 or 

2 



C 



d + 2 



(55) 



Then condition ( p4|) yields 



2D 

d< 2—D= d » 



(56) 



where d nc is the upper critical dimension in a melt [23]. The results (|55|)-(l56|) was 
obtained first in |E3|. 



Nevertheless, one can see that the exponent fl5"5| ) does not fulfill the self-consistency 
condition fl46| ) 



Q s t(x) oc J d D p 



1 — exp ( —i—px 



np\- 



(57) 



We will come to the contravention of the condition (f46| ) also in the case where ((d+2) > 
C(d + d° f + 2)-D. 



ii. The only way to satisfy the eqs. (f44D-(fl9|) is to impose on the exponents in eq. (|5T 
the condition 



C(d + 2)= C(d + d° + 2) - D > 2. 



(58) 



In this case (d°f — D = or 



C = Co 



2-D 



(59) 



On the other hand eq. (BSf) yields 
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d > 



2D 
2- D 



(60) 



i.e. the Gaussian solution ([59"D is self-consistent at d > d nc . In this case the entropic 
term dominates : {irp/N) 2 3> |Ei| oc [ | - The criterion (|60"D is equivalent to the 
embedding condition d°f < d and can be represented in the form 

2d 



D < D„ 



2 + d' 



The spectral critical dimension D s was discussed first in p3 



iii. At d — cL 



(61) 



Co(rf + 2) = Co(d + d° f + 2) - D = 2 



(62) 



and all terms have the same order : (irp/N) 2 oc |Ei| oc |S 2 |. The system is marginally 
stable and the stability condition is given by 

~4gL 



P 



- ci« c ,D) - c 2 (d uc D) 



> 



(63) 



which is always valid if pl d ^> 1. 



iv. If Qo{d + 2) = C{d + df + 2) — D < 2, i.e. d < d nc , the terns Si and £2 overwhelm 
the entropic one, (np/N) 2 |Ei| oc |E 2 |, and the system becomes unstable. Hence, 
at d < d nr the manifold is saturated in the melt, i.e. it loses its fractal nature and 



becomes compact [22 



B. Dynamic properties at d > d u 



We consider now the dynamics at d > d uc . There are two dynamic exponents, z and w. 
The exponent z measures the time dependence of a monomer displacement, i.e. 

ra+ioc fa r 



Q(t)= d D p 



2ni 



1-e 



St 



C(p, s) oc t 



2z 



(64) 
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and the exponent w measures the same for the center of mass 

C(p, s) oc t v 



ra+ioo d s 


1 - e st 


la-ioo 2%i 





(65) 



In eqs. (f>41)-(]5l)|) C(p,s) = (|R(p, s)| 2 ) is the Rouse-Laplace component of the correlator 
C(x,t) and the integral over s is taken along a straight line with Re(s) = a, so that the 
function is analytic at Re(s) > a. 



The formal solution of eq. (42) for C(p, s) is given by 

Cst(p) 



C(p, s) 



s + 



(66) 



£o+A/T(p, s ) 

where we have taken into account that at d > d uc the manifold is Gaussian, i.e. Q(p) = 
e(27ip) / N) 2 and C st (p) = l 2 (N/2np) 2 /Af. In eq. (|66| ) T(p, s) is the Rouse-Laplace component 
of the memory function fl34|). 

In order to calculate T(p, s) one needs analytically tractable approximations for the 
matrix density correlator Soo(k, t) in RPA and for the test manifold correlator -F(k; x, t). At 
{kl) d °fj\f < 1 we can use for the free system density correlator |16] 

FS\k,t) = FiSw[-k 2 D(k)t]\ (67) 

where D{k) = D /g(k), D = T/£ , g{k) is the Gaussian manifold static structure factor 
and is the corresponding static correlator. Using eq. ( pTD in the RPA-dynamic matrix 
(p3|), yields for the matrix density correlator 

S 00 (k, t) = S st (k) exp {-k 2 D coop {k)t} (68) 

where S Bt (k) is the static RPA-correlator, D coop (k) = p/£,oXst{k) and Xst'iM) = 



T 



F,t{k)\ + V(k). At (kl) i'N < 1 we have F£ ] w pM and Xst, 1 ~ V(k). Hence 

D coop (k) « f V{k). 



(69) 



,(0) 

At (fcZ) / A/" ^> 1 we probe only the local motion and, as a result, the dynamics is mainly de- 
termined by the single manifold behavior. So, the RPA-correlator Soo(k, t) is approximated 
by 



exp{-k 2 D coop (k)t} , (kl) d °fj\f < 1 
W = S 8t (k)x< Pi ; (70) 

where r = ^qI 2 /Td and z = (2 — -D)/4 is the Gaussian z-exponent. 
The corresponding Ansatz for the test manifold yields 

f expf-A; 2 ^*} , (kl) d °fj\f < 1 

where _Dg is the self-diffusion coefficient. 

The main problem now is to find the asymptotic behavior of T(p, s). If T(p, s) oc l/s 1_/3 
at s — > 0, where /? < 1, one should expect a renormalization of the Rouse dynamics. 

Let us consider first the dynamics in the intermediate displacement regime 

I 2 < Q(t) <R 2 G oc l 2 Af 2/d f. (72) 

We will show now that the asymptotic behavior of T(p, s) substantially depends on which 
limit in eqs. (|70|) -([7lD is relevant for the dynamics. 

i. Let us assume first that the dominant contribution for the integral over the wave 
vectors in eq. (0) comes from the interval 

AT 1 < (kl) d °f < 1. (73) 

Then we can use the second case in eqs. fl70|) and (71) as an input in the integral (p4|). 
By performing first the integration over wave vector k, and then over x, one can derive 
the result 

ATT® t) = /3|V(k = 0)| 2 ^^ ■ ' (74) 

where 

13 = z Q (d-d uc + 2). (75) 



19 



The Laplace transformation of eq. ([74]) at tqs <d 1 is given by 

AfT(P, s) = PW(k = 0)\ 2 Ssti t: 0) r^ flV"" . (76) 



l d + 

The condition f3 < 1 ( which is sufficient for the renormalized Rouse regime) immedi- 
ately defines the dynamical upper critical dimension 

AD 

due = = 2duc ( 77 ) 

i.e. the dimension above which the manifold has the simple Rouse behavior, at d = d uc 
we have the marginal Rouse behavior and only at d < d uc are the dynamic exponents 
z and w renormalized. The dynamical upper critical dimension d nc has been discussed 
first in [6,7], but the physical interpretation of this was different. It was asserted in 
ref. [6,7] that at d < 2d nc (or d < 2<ij) the strong entanglement effects become effective 
(see discussion in Appendix A). 

The substitution of eq. (|7|) in eq. ( j66|) and performing the inverse Laplace transfor- 
mation (by making use of the expansion theorem [|3I|) yields 

cfi ' H *^ WTi) (78) 

where A = T |^(k) | 2 5 s t(k = 0)/l d+2 and T(x) is the gamma-function. The eq. (|78|) 
is close to the stretched exponential form found by MC-simulation [ F28| . 

We should stress that the eq. flTBD was actually calculated in the limit p/N — > 0. That 
is the reason why we can use it first of all to comparison with MC-simulation results 
on the center of mass mean square displacement. By using eq. ([78]) in eq. (^) we 
obtain 

QUt) = ^-(-Y (79) 

where 



To 
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v ° - WTT) < 80 > 

and 

w = (5 = zo(d- d uc + 2). (81) 

We will compare the dependence ( fT5|) in the next Section with some recent MC findings. 
The using of eq. ([78]) in eq. yields 



z = z p = z 2 (d-d uc + 2). (82) 

Since (3 < 1, the condition z < holds, i.e. the interaction with the matrix slows 
down the monomer displacement. 



ii. If we assume now that the main contribution to the integral fl34j) comes from the small 

wave vectors 

{kl) d °sU < 1 (83) 

then by making use of the small wave vector approximation in eqs. (|70|) and (|7l|) we 
arrive at 

UY{p,s) (xs {d - d ^ )/2 . (84) 

Since <i > d uc , the simple Rouse behavior in the small wave vector regime does not 
change. 

Let us finally consider the large displacement regime 

Rl < Qcm(t). (85) 

In this case one should expect simple diffusive behavior 

Q cm (t) = dD G t (86) 
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where the self-diffusion coefficient ||f||| 



D G 



T 



(87) 



When the dynamics of the test manifold is characterized by the self-diffusion coefficient 
Dq, the dynamics of the matrix is driven mainly by the cooperative diffusion coefficient 
-Ccoop (see eqs. (ffO|)-([fl|)). Again, the small wave vectors interval, eq. (p3|), is relevant. 
Since in any case D coop S> Do, the calculations yields 



But D coop = O (N°) and d° f < d, therefore T(p = 0, s = 0) -> at Af -> oo. As a 
result Dq = T/A/"£o and again the simple Rouse behavior is not renormalized. 



in a melt composed of chains of the same nature. In order to calculate it, we have inte- 
grated over the melt's collective variables in the framework of the RPA and have used the 
self consistent Hartree approximation for the resulting effective action functional. It is very 
important that in this GRE the static and dynamic parts are treated in the same manner. 

In particular, if instead of the RPA the mode-coupling approximation (MCA) would be 
used for the collective variables (as it was done for linear polymer chains in 0,|J), one should 
substitute the RPA-correlators in eqs. (|34f) andfl3~7|) with the full correlators. It is then not 
obvious, how e.g. the simple screened form ( p8[) for the effective interaction potential could 
be obtained. In this respect generally speaking the dissipative and the reactive parts of the 
GRE in HIJ do not conform. 

The GRE derived here (as well as the GRE from [0-0]) cannot describe the reptational 
dynamics or entanglements. In order to do this, one should incorporate topological con- 
straints in the microscopic equation of motion. The extensive polymer mode-coupling ap- 
proach (PMCA) [4-7] which was worked out to explain the entanglement dynamics (without 




(88) 



IV. DISCUSSION AND CONCLUSIONS 



In the present work, we have derived the GRE (33)-(p7|) for the test polymeric manifold 
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the reptational model [16,17]) is, in our opinion, the result of misinterpretation of the GRE 
(see Appendix A). 

We have shown that at d > d nc the excluded volume interaction is screened out and the 
manifold is Gaussian with the exponent Co = (2 — D)/2 |22|,|23|]. Nevertheless the dynamic 



behavior is renormalized whenever d nc < d < d uc , where d uc , the dynamic upper critical 
dimension, is given by eq. (|77|), i.e. the reactive and the dissipative forces do not screen out 
simultaneously. 

For example, for the melt of polymer chains (D = 1) d nc = 2 and d uc = 4, and in the 3 - 
dimensional space one should expect the Gaussian static behavior but renormalized-Rousian 
dynamics. According to eqs. fl8i~D-(|82"D at d = 3 we obtain the exponents 



w = p = 3/4 = 0.75 (89) 

and 

2z = 3/8 = 0.375. (90) 

In eq. (|7^) we can take with a good approximation T(mf3 + 1) ~ T(m + 1), then the Rouse 
modes correlator has the stretched exponential form 

C(f,t) = CM exp 1-eA { (91) 



N J \r J 

The eqs. (91-93) were derived first in [4] as the solution of the first order approximation (so 
called RR- model). According to ref. [4-6] this regime describes the onset of the entanglen- 
ment dynamics (see Appendix A). 

As we have mentioned above eqs. (|74"D, (f?"6|) and (|78|) which lead to eq. ([H]) was actually 
calculated at p —>■ 0. That means that an experiment or simulation on the center of mass 
mean square displacement, Q cra {t), is the best candidate for the comparison with our results. 

MC-simulations of the bond fluctuation model [23J] as well as the MD-simulations p5|j2l 



of the athermal polymer melt have been undertaken. Recently also the static and dynamic 
properties of a realistic polyethylene melt have been studied by the extensive MD-simulations 
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27|,28|. Both in MC and MD simulations [24-28] a slowed-down motion at intermediate 



times for the center of mass mean square displacement is clearly observable. It was found 
e.g. that for the chain length N = 200 at the relatively short time (t < 3 • 10 6 MCS in 
Qcm(t) oc t w with w = 0.8 (instead of w = 1) in |24]] and w = 0.71 in P5] . This corresponds 
to our prediction w = 0.75. At larger times and scales the crossover to the reptational (for 
the long noncrossable chains) regime can be seen. This deviation from the simple Rouse 
regime also occurs for very short chains, N = 20, which clearly are not entangled. The 
regime Q C m{t) oc t ' 8 for the short-time regime has actually been observed first by Kremer 
and Grest pf. 

In Fig. 1 we have summarized the overall schematic behaviour for Q cm (t). At the rel- 
atively short times, r < t <C r R , and displacements, I 2 < Q cra {t) <C Rq, the test chain 
dynamics is mainly ruled by fluctuations from the interval (0), i.e. the Rouse dynamics 
is renormalized with the exponent (^). The picture which underlies this renormalization 
is visually represented in Fig. 2. The diffusing test chain in this case experiences mainly 
the short wavelength density fluctuations of the melt and, as a result, it is weakly "pinned" 
on the "lattice" induced by the density fluctuations. This "pinning" naturally results in a 
sub diffusive (w < 1) behavior at r < t <C r R . 

In the opposite limit, tr -C t and Rq -C Q C m(t), the long wavelength fluctuations from 
the interval ([83]) are relevant. Then the picture of the interplay between the test chain 
and the melt density fluctuations is given in Fig. 3. In this limit the melt almost does not 
influence the test chain and the simple Rouse regime is recovered. 

The crossover area, t pa tr and Q cm (t) pa R 2 1!) is not amenable for the theoretical in- 
vestigation mainly because of the lack of complete analytical expressions for the correlators 
(f?0|)-([7T|). Note that the renormalized curve in Fig. 1 converges asymptotically to the sim- 
ple Rouse curve from above. This is assured by the relationship T> > cITtq/C, , where the 
renormalized beads diffusion coefficient, Vq, is given by eq. (|80|). This condition can be 
represented in the form 
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In |24[] it has been explicitly shown that the renormalized beads diffusion coefficient (it is 



called an acceptance rate in lp4j| ) is larger than its Rouseian counterpart, i.e. the condition 
( |92|) holds for the real systems. The renormalized Rouse regime which is given in Fig. 1 is 
qualitatively the same as in MC [24] and MD [25-27] simulations (see e.g. Fig. 8b in [24], 
Fig. 9 in [27] and Fig. 3 in [28]). 

The slight stretching of the Rouse modes time correlation function has been found in 
[27,28] for modes with p — 1,2,3 (which still satiesfy the static p _2 -law). This is different 
from ref. [26] where no deviation from the ideal Rouse behavior for the higher modes has 
been observed. So as a future perspective it would be very interesting to solve eq. (44) 
numerically and compare the effects of the renormalization at finite p/N with the results of 
simulations. 

In [29,30] the MC-simulations for the dynamics of the athermal polymer melt have been 
undertaken. Special attention has beeb paid to the comparison of the crossable and non- 
crossable chains. These types of simulation are specially "designed" to check our results: 
The chains are long enough to follow the Rouse dynamics and to study its renormalization 
by interactions. It was shown in [29,30] that at relatively short wavelength Rouse modes, 
N/p < 6 and N = 500, the stretching parameter in the Rouse mode correlation function 
/? ~ 0.8. It is probably the result of the lattice structure structure influence. Unfortunately 
it is not clear why the static correlator of the Rouse modes deviates from the ideal val- 
ues even for the long wavelength modes (see e.g. Fig. 3 in [30]). For the relatively long 
wavelength Rouse modes, p/N > 10 and N = 500 (for crossable chains), the simple Rouse 
dynamics holds (i.e. (3 = 1 see Fig. 9 in [30]), but since the plot of Q cm {t) is not given 
explicitly, it stays inclear how the mode p — ► is renormalized. 

The temperature dependence of the renormalized Rouse regime is determined by the 
excluded volume interaction potential V(k). It is different from the static, where e.g. the 
screening effect does not depend on temperature at any potential. The condition for the 
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occurence of renormalization is given by £o ^ Io max N"T{p,t)dt, which, after using eq. (f74D , 
takes the form 



where t max oc « t N 2 . For the athermal case one can take as an estimate V(k = 0) « Tv 
(Edwards potential), and the condition (193|) always holds. For the pseudo-potential V ~ ev, 
where e has the dimension of molecular energy, the condition (j93|) is violated at rather high 
temperatures and the dynamics becomes Rousian. 
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Bundesministerium fiir Bildung und Forschung (BMBF) for financial support of the work. 

APPENDIX A: COMPARISON WITH SCHWEIZER'S POLYMER 



The purpose of this Appendix is to critically analyze the basic aspects of PMCA 0-|7j 
and to show that the pseudo-reptational exponents, which were obtained in PMCA, are a 
result of misinterpretation of the GRE. 

In Schweizer's PMCA the projection operator formalism and the mode-coupling approxi- 
mation were used in order to derive a GRE for the Gaussian chain. As a result the efffects of 
interaction of the test chain with the other chains are present in the form of a memory func- 
tion. The resulting GRE for the test chain correlation function C(a, (3; t) = (R(a, t)R(/3, 0)) 
takes the form Q 

2 We use here the notation of ref. Q 




(93) 
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f N f t Q Q2 

£ -C(a,ft;t) + J o d 1 J o dt'r(a,r,t-t')— 0(^,(3;?) = K S —C(a,ft;t) (Al) 
where K s = Zk^T ja 1 is the entropic spring constant and the memory function 



T(a,ft;t) = ( |U jT dk A; 4 (7 2 (k)c£; 2 (k)p m 5(k)F«(k, t) 



N 



x / d 1 d5Cj- l {a,ft;\^Cj Q {^,8;^,t)Cj' 1 {8,ft;\i) (A2) 



o 



where p m is the average density, C*(k) is the direct correlation function, u)(k) and ^(k) are 
static correlators for the test chain and matrix correpondingly, uj®(a, ft; k, t) is the dynamic 
test chain density correlator associated with segments a and ft, oj(a, ft; k) = u Q (a, ft; k, t = 
0) and iv(k, t) is the normalized dynamical collective matrix correlator. As usual the 
superscript Q denotes the evolution via projected dynamics. 

As a main approximation the projected single chain dynamics was substituted by so- 
called " renormalized Rouse" (RR) dynamics 



uj Q (a, ft; k, t) -> u RR (a, ft; k, t) w w(ot, /5; ]t)F RR (k, t). (A3) 



This RR-dynamics is a first order approximation of an iterative solution of eq. (|A1|) when 
as a zeroth order approximation the well-known Rouse (R) expression 



cu Q (a, ft; k, t) -> w fl (a, /3; k, t) w w(a, /3; k)F R (A;, t) (A4) 



is employed. The corresponding justification given in 0-0] to take the i?i?-approximation 
as a starting point "in crude analogy with Enskog theory" remains questionable. 

For the projected collective correlation the real (not RR) dynamic evolution is employed 

F?{k,t)^F s {k,t) (A5) 



This approximation could be justified for small k [|38| and was used e.g. in the glass transition 
theory |39 . 



The main problem with the present analysis of eq. ( |A1| ) is that it cannot be solved 
iteratively. It is easy to see that eq. (|A~I) is substantially non-linear because in the memory 



function (A2) the dynamic correlator is given by 
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u Q (a, (3; k, t) = (exp {-ik [R(a, t) - R(/3, 0)]}) 

= exp | -y [C(a, a; t = 0) - C(a, (3; t)] j . (A6) 

The first line in eq. (|A£|) implies that the real dynamical evolution is used and the second 
equality comes from the fact that the fluctuations of a .R-variable are Gaussian. As a 
result eq. (|A1|) is substantially non-linear and should be treated self-consistently around the 
bifurcation points as in the glass transition theory |39| . Substantially non-linear means that 
the range of parameters (temperatures, length of chains, time, etc.) is such that the mode- 
coupling friction term in eq. (|A1| ) is comparable to or much larger than the bare friction 
one. But this is exactly the range which PMCA investigates as a crossover from Rouse to 
entangled dynamics. 

Instead of selfconsistent solutions around a bifurcation point iterations are actually used. 
As a zeroth order approximation this yields the .R-expression, and the first order approxima- 
tion - the RR-model - is used as a starting point for the description of entanglements. The 
second order approximation gives already the result which looks like entangled dynamics. 
But one can see that this iteration procedure is divergent as it should be for a substantially 
non-linear equation. 

For example, the Rouse friction coefficient £r oc D^ 1 oc A, the -RR-model leads to 
£rr oc D~^ r oc A 3 / 2 and at last the next approximation yields £ coi i oc D~ oil oc A 2 . 

There is a number of differences between Schweizer's GRE and the GRE which was 
derived here. First of all, the statics and dynamics of the matrix (or melt) in eqs. (35)- 
(41) are described by RPA-correlators. This results in the conventional static screening 
(see eq. (40)). In order to introduce full correlators we should go beyond the Gaussian 
approximation in the expansion (17), i.e. take into account the melt fluctuations around the 
RPA [34]. This is a tough problem even for statics and it is not clear beforehand whether 
the simple bilinear form of the memory function (36) will remain. Unlike this, in Schweizer's 
GRE [4] the full melt correlator appears in the memory function whereas for the static part 
it is simply apriori taken to be the screened form of the interaction (see eq. (2.3) in ref. [4]). 
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Secondly, by the analytical investigation of Schweizer's GRE different analytically 
tractable formulas for the matrix correlator were assumed. For example in sec. Ill c of 
ref. [4] ("Center of mass diffusion and shear viscosity") or in sec. Ill in ref. [5] the frozen 
matrix i.e. F s (k,t) = 1 was adopted. This brings the result for the RR-model Dq oc N~ 3 ^ 2 , 
which was mentioned above. In our GRE we cannot consider this case because in RPA- 
approximation the matrix is not frozen. But in most of the computer simulations [24-29] 
the athermic polymer melts have been studied, which cannot be frozen. In this case at small 
wave vectors ({kl) d °fj\f <C 1) the dynamics of the matrix is driven by -D COO p and the simple 
Rouse dynamics of the test chain is recovered. 

For the small or intermediate displacements, I 2 < Q cm (t) < R%, the so called Vineyard 
approximation, i.e. F s (k, t) = F(h, t), was adopted. This leads for the RR-model (first order 
approximation) to formally the same mathematics as in our sec. Ill B(i). As a result the 
values of the exponents w and z (see eqs. (91)-(92)) as well as the dynamical uppper critical 
dimension, d nc = 2d°p are the same. This value of d nc is assured by the short range nature of 
the interaction function V(r) in the memory kernel which leads to the scaling for the number 
of dynamically effective contacts: Z^yn oc N 2 ~ d '/ d °f . This scaling law has been discussed first 
in ref. [6,7]. Nevertheless the criterion d < 2d°j has been considered as a necessary condition 
for the onset of the entangled dynamics. As we have already mentioned above in order 
to obtain the entangled dynamics exponents in PMCA the next (second order) iteration is 
implemented. In the spirit of PMCA we could use e.g. the renormalized Rouse exponents 
as an input in the r.h.s. of eqs. (83)-(84) and rederive the pseudo-reptational exponents 



W- 



rep 



rep 



z{d - d uc + 2) 
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(A7) 
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These exponents have been given in |J. Unfortunately, this way of analysis does not seem 



to bring reliable results. 
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In the present work we have made every effort to prove that the excluded volume inter- 
action results in the renormalized Rouse regime in the melt. We believe that the topological 
constraints are essential for the entangled dynamics. The recent MC-simulation p9| , |30f have 
shown the decisive role of the topological constraints underlying the reptational dynamics. 
In particular the crossable and noncrossable models of rather long chains (10 < N < 500) 
were used for the simulation of statics as well as dynamical properties of polymer melts. It 
was shown that the static properties of both models are absolutely identical (see e.g. Fig. 3 in 
||29||). On the other hand the dynamic behavior is completely different. The crossable chains 
(irrespective of short or long) show at the relatively long wavelengths the simple Rouse be- 
havior (see e.g. Fig. 6 and Fig. 9 in |30]]). Also the self-diffusion coefficient has the Rouseian 
scaling Dq oc N~ l (Fig. (7) in p9|). For the noncrossable chains (at N = 100,300 and 
500) the stretched exponential regime is found at all mode wavelengths. The self-diffusion 
coefficient scales as Dq oc iV~ 2 ' 08 at iV > 40 and the mean-square displacement Q(t) oc t 1//4 
in the corresponding time window. 

This proves that the chain-crossing condition does not touch the static properties but 
has a dramatic effect on the dynamics. On the other hand, the static input completely 
predestines the dynamic behavior in the PMCA-formalism. This contradiction unfortunately 
put severe doubt on the PMCA. We rather feel that the explicit taking into account of 
topological constraints in the microscopical equations of motions is absolutely required. 
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FIGURES 
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FIG. 1. A schematic plot of Q CYa {t) for the simple Rouse (dashed line) and the renormalized 
Rouse (solid line) dynamics 
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FIG. 2. The test chain experiences the short wavelength melt density fluctuations and as a 
result it is weakly "pinned". 
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FIG. 3. The interplay of the test chain and very long wavelength density fluctuations does not 
influence the Rouseian dynamics. 
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